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Abstract. We present two models for the cosmological UV background light, and calculate the opacity of GeV gamma-rays 
out to redshift 9. The contributors to the background include 2 possible quasar emissivities, and output from star-forming 
galaxies as determined by recent a semi-analytic model (SAM) of structure formation. The SAM used in this work is based 
upon a hierarchical build-up of structure in a ACDM universe and is highly successful in reproducing a variety of observational 
parameters. Above 1 Rydberg energy, ionizing radiation is subject to reprocessing by the IGM, which we treat using our 
radiative transfer code, CUBA. The two models for quasar emissivity differing above z = 2.3 are chosen to match the ionization 
rates observed using flux decrement analysis and the higher values of the line-of-sight proximity effect. We also investigate 
the possibility of a flat star formation rate density at z > 5. We conclude that observations of gamma-rays from 10 to 100 GeV 
by Fermi (GLAST) and the next generation of ground based experiments should confirm a strongly evolving opacity from 
1 < z < 4. Observation of attenuation in the spectra of gamma-ray bursts at higher redshift could constrain emission of UV 
radiation at these early times, either from a flat or increasing star-formation density or an unobserved population of sources. 
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1. INTRODUCTION 

The Extragalactic Background Light (EBL) is the to- 
tality of light emitted by galaxies over the history of 
the universe. High-energy gamma-rays propagating over 
cosmological distances are attenuated due to electron- 
positron pair production with EBL photons. The optical 
depth tends to increase with gamma-ray energy, which 
modifies the high-energy spectra of distant sources. 
The cross section for this annihilation process peaks 
at a wavelength of 1.24(E r /TeV) microns. For Fermi 
(GLAST) energies (< 300 GeV), the relevant wave- 
lengths are in the ultraviolet, and for energies below 
about 100 GeV the ionizing background is responsible 
for the pair-creation interaction. A previous calculation 
of the background in the UV found that the contribution 
from resolved galaxies could make the universe optically 
thick to —30 GeV photons [1]. 

In the optical and near-UV, the EBL is constrained by 
firm lower limits from number counts, e.g. [2, 3]. Indi- 
rect upper limits have recently been placed on the optical 
and near-IR background from ground-based gamma-ray 
experiments e.g. [4, 5], However, little is known about 
the background at energies above the Lyman limit from 
direct measurements. Most ionizing photons from star- 
forming galaxies are absorbed by local cold gas and dust, 
with an uncertain fraction f esc escaping to the intergalac- 
tic medium. Predicting optical depths to lower energy 
gamma-rays is further complicated by the fact that the at- 
tenuation edge for an optical depth of unity increases to 



a redshift of several, meaning that the evolution of ion- 
izing sources must be understood to high redshift. Be- 
yond a redshift of about 2, the faint-end of the quasar 
luminosity function remains highly uncertain, with lim- 
its placed on the ionizing background via observations of 
the hydrogen and helium Lyman opacities. The rapid in- 
crease of neutral helium above z— 3 is believed to signal a 
strong decrease in the quasar emission of hard spectrum 
UV radiation. While stars are thought to be the domi- 
nate source of UV photons at higher redshift, uncertainty 
in star-formation rates and efficiency, and possibility of 
changing escape fraction or initial mass function make 
predicting optical depths for gamma rays at high redshift 
much more difficult than studies of absorption with local 
TeV sources. 



2. MODELING 

We have calculated the EBL from the infrared to soft X- 
rays out to redshift 9.5 for 3 scenarios varying in escape 
fraction, quasar emissivity, and high redshift star for- 
mation rate. Our approach utilizes recent semi-analytic 
modeling work [8] to determine stellar emissivities from 
galaxies. These models are based on ACDM with h = 
0.7, £l m = 0.3, and use prescriptions based upon N-body 
and hydrodynamic simulations together with normaliza- 
tion from low-redshift data. The models contain all the 
ingredients thought to be important in the evolution of 
galaxies. A successful model of the UV background must 
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FIGURE 1. The background (from x-rays to IR) for the fidu- 
cial (black) and quasar-dominated (gray; green in online edi- 
tion) models at the indicated redshifts. The backgrounds at 
nonzero redshift are for photon populations evolved to present- 
day. The data include lower limits from GALEX, HST, and 
IRAC (upward pointing arrows). Direct detection data shown 
with other symbols include HST and DIRBE data with subtrac- 
tion of zodiacal light and galactic sources. 
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FIGURE 2. The total emissivity from star-forming galaxies 
in our semi-analytic model [8], at the given redshifts, scaled 
to 1 MQyr~ l . HI and Hell ionization energies are indicated by 
thin vertical lines. 



take into account reprocessing of ionizing photons by the 
intergalactic medium (IGM). Neutral hydrogen and neu- 
tral and singly -ionized helium in IGM clouds act as both 
sinks and sources of ionizing radiation [9] . We have used 
CUBA [10] (an update to this software is in preparation), 
a numerical code for determining the propagation of UV 
photons through the IGM, to calculate the background 
as well as relevant observable parameters such as ioniza- 
tion rates and fractions. The quasar luminosity evolution 
is based upon the 3 models presented in [7] (see Figure 
4). We assume a power law quasar spectrum with slope 
1.57 shortward of 1200 Angstroms, as suggested by HST 
studies of quasarsfll]. We have developed 3 models for 
the evolution of the UV background that explore possi- 
bilities in star formation and quasar output; see Figures 1 
through 4. 

For our fiducial model, we have used the 
observationally-based quasar evolution model of [6], 
with a moderate f esc of 0. 1 . It should be pointed out that 
this model produces emissivities that are quite similar to 
those of model 'A' in [7]. In this case, stars and quasars 
produce enough ionizing photons to match ionization 
rate data from analysis of Lyman alpha flux decrement. 

The quasar-dominated model is based upon quasar 
model 'C in [7], and is motivated by the higher ion- 
ization rates inferred from measurements of the line-of- 
sight proximity effect. We also use a minimal f esc of 0.02 
in this case. 

A third possibility we have considered in our fiducial 
model is to modify the steeply declining comoving star- 
formation rate density in our semi-analytic model to 
remain flat for z>5. 



3. DISCUSSION AND CONCLUSIONS 

We have presented models of the UV background at 
all wavelengths relevant to absorption of gamma-rays in 
the Fermi LAT energy range (0.02 to 300 GeV). Ma- 
jor uncertainties in UV emission include the poorly con- 
strained escape fraction from star-forming galaxies and 
the faint end of the quasar luminosity function above 
z ~ 2, and we have explored the resulting EBL and and 
consequent gamma-ray absorption out to redshift 9 in 
three cosmological scenarios varying in these parame- 
ters. Our fiducial model is based upon the quasar lumi- 
nosity functions of [6]. While the faint end of the quasar 
luminosity function remains unobservable at these red- 
shifts, it is argued in this work that bright quasars at the 
primary contributor to the total quasar emissivity in this 
epoch, and the contribution from unseen faint quasars 
is small. The ionizing flux from quasars must also de- 
crease rapidly above z^3 to account for sharply increas- 
ing helium Lyman-opacity, and must constitute only a 
small fraction of the total ionizing flux by z^6 [12]. 
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FIGURE 3. The cumulative emissivity density of quasars at 
912 (upper) and 228 (lower) Angstroms in our models. The 
black curve is our fiducial model, based on [6], and the dashed 
gray/green curve is based upon the model "C in [7] multiplied 
by a factor 0.8. 



FIGURE 5. The ratio of column densities Hell/HI (softness) 
arising in each of our models. The much more rapid rise for our 
fiducial model is indicative of its faster transition to a stellar- 
dominated ionizing background above z~2.5. 
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FIGURE 4. The star formation rate density in our semi- 
analytic model. The dashed line showes the flat SFR density 
we consider as a possible modification to the declining SFR 
density of our model [8]. Data is scaled to our Chabrier initial 
mass function. 



With a moderate contribution to the background from es- 
caping radiation from star-forming galaxies, this model 
produces good agreement with UV softness measure- 
ments (Figure 5) and the low ionization rate (Figure 
6) inferred from recent flux decrement measurements, 
cf. [13]. Our 'quasar-dominated' model is based upon 



model 'C in [7] and leads to the more absorption at low 
and intermediate redshifts. Inclusion of this scenario is 
motivated by the ionization rate calculated from obser- 
vations of the line-of-sight proximity effect, which are 
typically 2-3 times higher than those derived by other 
methods. A large stellar ionizing component at this flux 
would violate observational limits on f esc [14]. A num- 
ber of potential biases exist in this calculation which 
could account for the the higher r_i2, though the flux 
decrement method is not without biases either [15]. The 
quasar-dominated scenario also has trouble reproducing 
the trend of increasing UV softness approaching the red- 
shift of helium reionization. 

As an additional possibility, we have modified our 
fiducial model to have a flat star formation rate at z > 5. 
Predictions for gamma-ray absorption in these 3 scenar- 
ios are shown in Figure 7. All of our models create an 
optical depth that changes greatly between z = 1 and 
4. The modified SFR density at z > 5 leads to signifi- 
cantly greater optical depth at these redshifts, an effect 
that may be discernable if sufficient high-z VHE sources 
can be observed. The mean free path of ionizing radiation 
is expected to fall quickly above z^6, as suggested by 
the appearence of complete Lyman alpha Gunn-Peterson 
troughs [16]. However, the (1+z) multiplier from ob- 
served to rest-frame energy means that photons observed 
in the 10-100 GeV energy decade are affected primar- 
ily by non-ionizing UV and optical photons at these 
high redshifts. This means that significant absorption of 
gamma-rays in this decade of interest can occur due to 
radiation below the Lyman limit even during reionization 




FIGURE 6. The ionization rate T_i2 of hydrogen in units 
of 10~ 12 s , for each of our models. As before, our fiducial 
model is in black, the quasar dominated model in dashed gray 
(green). Data in black is from flux decrement measurements of 
the Lyman alpha forest. The generally higher gray/green points 
are from measurements of the line-of-sight proximity effect, 
and motivated our consideration of the higher emissivities in 
the second model. 



FIGURE 7. The attenuation edges for gamma-rays in each 
of our models. The fiducial and quasar-dominated models are 
indicated using solid and dashed lines respectively, and the 
dotted lines show the model with a flat SFR density above z = 5 
(Figure 4). Each set of curves indicates the redshift at which the 
optical depth T reaches the indicated value (0.5, 1, 2, 3, and 5, 
from lower left to upper right). 



when the universe is largely opaque to ionizing UV pho- 
tons. While most data from the highest redshift Lyman- 
break galaxies has suggested a declining SFR [17], other 
techniques such as measurement of gamma-ray burst dis- 
tributions [18] have found higher rates, and little data ex- 
ists at these redshifts. An upcoming publication will ex- 
amine different high-redshift scenarios in much greater 
detail, and in the context of reionization. 
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